Order Parameter to Characterize Valence-Bond-Solid States in Quantum Spin Chains 



O 

o 

(N 
o 

Q 

o 



o 



CO 

o3 



i 

c 

o 
o 



> 



O 



-a 
c 

o 
o 



X 
S3 



Masaaki Nakamura 1 and Synge Todo 2 ^ 3 
1 Department of Applied Physics, Faculty of Science, 
Science University of Tokyo, Kagurazaka, Shinjuku-ku, Tokyo 162-8601, Japan 
2 Theoretische Physik, Eidgendssische Technische Hochschule, CH-8093 Zurich, Switzerland 
3 Institute for Solid State Physics, University of Tokyo, Kashiwa 277-8581, Japan 

(Dated: February 1, 2008) 

We propose an order parameter to characterize valence- bond- solid (VBS) states in quantum spin 
chains, given by the ground-state expectation value of a unitary operator appearing in the Lieb- 
Schultz-Mattis argument. We show that the order parameter changes the sign according to the 
number of valence bonds (broken valence bonds) at the boundary for periodic (open) systems. This 
allows us to determine the phase transition point in between different VBS states. We demonstrate 
this theory in the successive dimerization transitions of the bond-alternating Heisenberg chains, 
using the quantum Monte Carlo method. 

PACS numbers: 75.10. Jm,75.40.Cx 



About two decades ago, Haldane made a striking con- 
jecture that the integer-spin 5 antiferromagnetic Heisen- 
berg chain has a unique disordered ground state with an 
energy gap, while for a half-odd-integer spin there is no 
energy gap and the system belongs to the same univer- 
sality class as the 5=1/2 case For 5 = 1, this con- 
jecture was confirmed by many numerical |^] and exper- 
imental studies || . Affleck, Kennedy, Lieb, and Tasaki 
studied an 5 = 1 isotropic spin chain with special bi- 
quadratic interactions which has the exact ground state 
with a finite gap Q|. They also proposed the valence- 
bond-solid (VBS) state for the Haldane gap systems, and 
concluded that the ground state of the Heisenberg chain 
is described approximately by the VBS state. Since the 
spin configurations of the VBS state show the hidden an- 
tiferromagnetic order, den Nijs and Rommelse proposed 
the string order parameter to characterize the 5=1 Hal- 
dane phase ||, 
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where a = x,y, z and l^o) means the ground state. Thus 
this order parameter enables us to detect the VBS state 
indirectly. The string order parameter was generalized to 
5 > 1 cases by Oshikawa (6). 

On the other hand, Haldane's prediction was also dis- 
cussed by the Lieb-Schultz-Mattis (LSM) type argument 
0. Using this technique, Affleck and Lieb examined Hal- 
dane's conjecture |@). Later, Oshikawa, Yamanaka, and 
Affleck generalized this argument, and obtained a nec- 
essary condition for a gapped state Q. However, the 
relation between the LSM argument and the VBS pic- 
ture including the string order parameter has not been 
fully understood. In this Letter, we discuss this relation, 
and show that an overlap integral appearing in the LSM 
argument [Eq. (||) below] plays a role of an order pa- 
rameter which detects VBS ground states directly. We 



also show that it can also be applied to determine phase 
boundaries in between different VBS states. We demon- 
strate this idea by the quantum Monte Carlo (QMC) 
simulation for successive dimerization transitions in the 
bond-alternating Heisenberg spin chains. 

First, let us review the LSM argument briefly based on 
Ref. [[| and introduce the order parameter. We consider 
a periodic spin chain of length L with short range inter- 
actions. We assume the translational (T : S" — ► 5" +1 ) 
and the parity (V : 5" — > 5g +1 _.-) invariance. Now we 
introduce the following "twist operator" 



U = exp 



i=i 



(2) 



Since this operator rotates all the spins about the z- 
axis with relative rotation angle between the neighboring 
spins 2n/L, it creates spin- wave- like excitations. Apply- 
ing the twist operator q times to the unique normalized 
ground state |^o) generates a set of low- lying excited 
states {|*fc) = ?7 fc |*o)} (k = 1, • • ■ , g). The excitation 
energy of the state is evaluated as AE ~ C(L _1 ). 
Although is not necessarily an eigenstate of the 

Hamiltonian, there exists at least one eigenstate with en- 
ergy of C(L _1 ), if the state l^i) is orthogonal to the 
ground state |^o)- 

In order to consider the orthogonality of these states, 
we introduce the following overlap integral which plays a 
central role in this Letter: 

4 9) ^ <*o|* g > = <*o|^|*o>- (3) 

The invariance under transformations T and V yields 

4 9) = {^^TVT-^o) =tP*( s i-™)*z$\ (4) 
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where m is the magnetization per site. Eq. fl) shows 



that 



7^ is possible only when 5 — m = integer. 



2 



In this case, the system has a gap without breaking the 
translational symmetry. On the other hand, for a rational 
value S — m = p/q with p being an integer, |Wq), l^i),- • •, 
l^g-l) are mutually orthogonal. This means that the sys- 
tem in the L — ► oo limit is gapless {z£} = 0), otherwise 
gapped with q-fold degenerate ground state (z^ 7^ 0) 
due to the spontaneous breaking of the translational sym- 
metry. Thus q(S — to) = integer is a necessary condition 
for a gapped state [_ 
imaginary due to Eq. (||) 

In the previous works based on the LSM argument 
|8[ ||, however, explicit values of Zl have not been 
calculated. Now let us evaluate zl for various ground 
states. In the Neel state with two-fold degeneracy, one 
can obtain zf^ = (— 1) 2S immediately. In order to calcu- 
late Zl in VBS states, we introduce the Schwinger boson 
representation for the spin operators |l0|: 



In this case, zi is real or pure 
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where these bosons satisfy the commutation relation 
[ai,aj] = = £>ij with all the other commutations 

vanishing. The operator (&j) increases the number 
of up (down) S — 1/2 variables under symmetrization. 
Then, a generalized VBS state discussed in Ref. in a 
periodic system (see Fig. 0(a)) is written as 
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where B] , = alb] — b]al, J\f is a normalization factor, 
and |vac) is the vacuum with respect to bosons. The 
integers to and n satisfy to + n = 2S. Using relations 
Ua^U' 1 = a]e+' l ^/ L and Ub)U- x = fe]e~ i7r ^ L , a twisted 

valence bond UBj J+1 C/ — 1 for 1 < j < L — 1, and that 
located at the boundary are calculated as follows, 

'•'CV ' - e-Wi J 6 j +1 _^ 6 J J +ij (8) 
UB^U- 1 = -(e-^ L a L b\-^ L b{a\). (9) 

In the latter case, a negative sign appears for each valence 
bond. This reflects a property of an S = 1/2 spin which 
needs 4ir rotation to return to the original state. Thus 
the asymptotic form of zl is given by 



(m,n) I 
VBS 



\U\^) = (-ir[l-0(l/L)}. (10) 



It turns out that zl changes its sign according to the 
number of valence bonds at the boundary. The value 
of zl for to = n — 1 was calculated by Totsuka and 
Suzuki |l^] . A similar relation was also found by Bones- 
teel who discussed two-dimensional S — 1/2 dimer sys- 
tems jD^. In Eq. (|l0|), the factor (—1)" originates form 
the relative twist angle 2ir at the boundary, so that it 



(a) 



(b) 



FIG. 1: Generalized VBS states with (m,n) = (4,2) for (a) a 
periodic chain and for (b) an open chain with even L. Broken 
ovals mean the symmetrization of S — 1/2 variables at each 
site, and solid triangles denote isolated bosons. 



does not depend on the detailed structure of the twist 
operator (||). Therefore, the definition of zl is still mean- 
ingful even in cases with broken translational symmetry 
(to 7^ n) which does not accord with the LSM argument. 
Especially, in the present definition, zl has a useful sym- 
metry zl — * (— 1) 2S zl under interchange of to and n. 

In order to calculate Zl in more detail, we introduce 
the matrix product (MP) formalism Jl4|, [hq] . This is 
useful to calculate the ground-state expectation values 
and correlation functions. In this formalism, Eq. (0) is 
expressed as 
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These matrices for general integer S case with to = n are 
given in Ref. jl2|. Based on this result, the matrices for 

the twisted VBS state C/'I^vbs^ are S iven b y 
gqj(s,r) = (_i)«-s+V("-2s+2)W£ 



x v / m C r _ ln C s _ 1 ( a t) m - r+s (6t)"+'- a |vac) J -, (12) 



x VmC r _i n C s _i(at) m - r+s (6])"+'- a |vac> J -, (13) 

where l<r<m + l, 1 < s < n + 1, and |vac).,- is the 
boson vacuum at the j-th site. By introducing a transfer 
matrix defined as 

Gf(r 1 ,r 2 ;s 1 ,s 2 ) = 9 j(ri,s 1 )gg J (r 2 ,s 2 ), (14) 

we obtain zl as 



= (-1)'' 



t Tr[G^fF 2 
Tr[G£G$]W 



(15) 



Although the analytic form of zz, is complicated in gen- 
eral, one can evaluate numerical values of Eq. (|l5|), which 
will be shown in Fig. below. In the fully dimerized 
case (m,n) = (2S 1 , 0), zl is given by a simple form 



v (i) 



Efcto ei(2 ^ 2fe)7r/L /(2^+l)] L/2 - 



One can show that zl behaves similarly even in open 
systems. For even L (see Fig. 0(b)), the VBS state of 
Eq. (|) has (n + l) 2 -fold de generacy due to the isolated 
bosons at both ends Q. Then, zl changes its sign ac- 
cording to the number of isolated bosons at the L-th site, 



(-l) n e- lror / z [l-0(l/L)]. 



(16) 
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FIG. 2 

with S = |, 1, |, and 2, obtained by the QMC calculation 
The horizontal lines indicate zl for the VBS states calculated 
by the MP method [Eq. @]. 



of the L = 64 periodic BAHC FIG. 3: System size dependence of the successive dimerization 



transition points <5 C (L) of the periodic BAHC with S = |, 1, 
|, and 2. The extrapolation to the L — > oo limit is done by 
5 C (X) = 5c(oo) + A/L 2 + B/i 4 + C/L G (solid lines). 



Although is complex in this case, the imaginary part 
vanishes in the large-!/ limit, and the real part of zl 
behaves in the same way as that of the periodic cases. 
For open systems with odd L, the number of isolated 
bosons at the L-th site is to, and the ground state has 
(to + l)(n + lVfold degeneracy. Therefore, zl is simply 
given by Eq. ( |l6| ) multiplied by (— 1) 2S . 

Next we consider zl in connection with a low energy 
effective theory. According to Schulz's bosonization anal- 
ysis Jig , the Lagrangian density of the uniform Heisen- 
berg spin chain is given by the sine-Gordon model, 

C= ^K [ V< ^' T )] 2 - ^™s[qV2cj>(x,T)}, (17) 

where r is the imaginary time, a is a short range cut 
off, and K and y<f, are the parameters determined phc- 
nomenologically. The phase field is related to that used in 
Ref. @ by = 2VsV>i, and q = 1 (q = 2) for S integer 
(S half-odd integer), where q corresponds to degeneracy 
of gapped ground states. In the gapped (gapless) region 
one has y${l) — > ±oo (j/<^(0 — > 0) for I — > oo under renor- 
malization a — > e'a. On the unstable Gaussian fixed line 
[y<p(0) = OwithXq 2 < 4], a second-order "Gaussian tran- 
sition" takes place between the two gapped states. In this 
formalism, and the spin wave excitation created by U cor- 
responds to the vertex operator exp(iv / 20), so that zl for 
q(S — to) = integer is related to the ground-state expec- 
tation value of the nonlinear term as zf^ oc (cos(g-\/2</>)) 
and the three fixed points — ±oo,0 correspond to 
Zoo = Tl,0, respectively. Thus the Gaussian critical 
point can be identified by observing = p7| . 

In order to demonstrate the above argument, we con- 
sider the bond-alternating Hcisenberg chain (BAHC), 

L 

n = jJ2^-s(-iy}s r s J+1 . (is) 

i=i 



For this model, the VBS picture is considered to be re- 
alized approximately: The configuration of the valence 
bonds (to, n) changes from (0, 2S) to (25, 0) successively 
as 6 is increased form —1 to 1, meaning the existence of 
2S quantum phase transitions ||, [ij] . Since the transla- 
tional symmetry is explicitly broken, the effective model 
of the BAHC is given by Eq. (0) with q = 1 for any S 
pUf , and these transitions are described as a Gaussian 
type. 

In Fig. |[ we show z^ (with q = 1) of the L = 64 peri- 
odic BAHC with S = 1, |, and 2 as a function of 5. For 
the calculation of zl, we employ the QMC method with 
the continuous-time loop algorithm ]l9| . We used the 
multi-cluster variant of the loop algorithm. The QMC 
steps are 10 3 for thermalization, and 5 x 10 5 -10 6 for the 
measurement, and the inverse temperature (3 is taken 
large enough so that the value of zl can be identified as 
that of the ground state. The largest inverse temperature 
used in this calculation is (3J = 128. 

As seen clearly in Fig. ^, the 5 dependence of z^ agrees 
qualitatively with the present interpretations based on 
the VBS picture, which predicts z^ — (— l) n with a sym- 
metry zl — > (— 1) 25 zl for 5 <-> —5 (to <-» n). We have 
also calculated zl for open chains with even and odd 
L's, and confirmed that the results agree with our predic- 
tions including Eq. Jl6|). In Fig. |^, we also present zl for 
the VBS states calculated by the MP method [Eq. @] 
as horizontal lines. The difference between zl for the 
BAHC and that for the VBS states becomes larger as 
S increases. This means that the VBS picture becomes 
poor for the BAHC with large S. 

Next, we determine the critical point i5 c by observ- 
ing zl = with L up to 320 (Fig. |^). In this calcu- 
lation, the inverse temperature is taken as (3 J = L/2S, 
being assumed the Lorentz invariance. Extrapolation to 
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the L — > oo limit has been done by the least-squares 
fitting by assuming the asymptotic form as 8 C (L) = 
S c (oo) + A/L 2 + B/L A + C/L 6 . This polynomial with 
even powers of 1/L is justified by the parity symmetry 
which ensures that Zl is real [Eq. (||)]. For S = 1, we 
have obtained the critical value for the transition be- 
tween the (1,1) and (2,0) phases as S c = 0.25997(3), 
where ( ) denotes 2a. This result is consistent with the 
previous estimates: S c — 0.2595(5) by the QMC calcu- 
lation for the susceptibility |2(J and 0.2598 by the level- 
crossing method ]2l[ |. Similarly, we identify the criti- 
cal point of the S = 3/2 case [(2, l)-(3, 0)] and those 
of the S = 2 case [(2, 2)-(3, 1), (3, l)-(4, 0)]. We obtain 
S c = 0.43131(7), 0.1866(7), and 0.5500(1), respectively 
(see the inset of Fig. ||). They are also consistent with 
5 C = 0.4315, 0.1830, and 0.5505 obtained by the level- 
crossing method [p2| , but much more accurate. 

Here, we comment on the method used in Refs. plj] 
and p2j proposed by Kitazawa who pointed out that 
the Gaussian transition with q = 1 can be identified 
by a level crossing of excitation spectra under twisted 
boundary conditions (2^]. Since this method is also ex- 
plained by the sine-Gordon theory |2i| and the VBS pic- 
ture [^l[ |22j , the results are considered to be equivalent to 
ours. Although finite-size corrections in the level-crossing 
point in VBS states tend to be smaller than those of the 
Zl = point, the application of the level-crossing method 
to larger systems is difficult, because it relies on the exact 
diagonalization to obtain excitation spectra which needs 
the whole Hilbert space. On the other hand, our ap- 
proach is based only on the ground state quantity, Zl, 
so that various numerical methods such as the present 
QMC method can be employed. This makes it possible 
to deal with systems with enormous Hilbert space, such 
as large-.? and ladder systems. Note that the density ma- 
trix renormalization group method is also suitable for our 
approach, since zl remains as a meaningful order param- 
eter even in open systems. Furthermore, our method also 
has an advantage to comprehend the physical picture as 
already discussed. 

Finally, we discuss the relation between zl and the 
string order parameter. In Ref. ||, Eq. (Q) for I^vbs^) 
was calculated as 0* tring = [ " ( +"^ ] 2 <5n,odd, where finite 
and vanishing Of tring correspond to negative and positive 
Zl, respectively. However, this calculation is limited only 
for S integer, and the twist angle 7r in Eq. ([j]) is needed 
to be generalized for other cases [||. In addition, it is 
difficult to determine accurate phase transition points by 
O. 
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because it does not change the sign. Thus zl 
turns out to be more rational order parameter to describe 
VBS ground states in a unified way. 

In summary, we have introduced zl given as the 
ground-state expectation value of the twist operator (||) 
as an order parameter to characterize various ground 
states in quantum spin chains. Especially, zl changes 



the sign according to the configuration of valence bonds. 
This property enables us to determine the critical point 
between different VBS states by observing zl = 0. We 
have demonstrated this theory by using the QMC sim- 
ulation for the successive dimerization transitions of the 
BAHC, and determined the phase boundaries with quite 
high accuracy. 
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